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DIFFUSION OF HEAT FROM A LINE SOURCE IN ISOTROPIC TURBULENCE! 


By MAHINDER 8. Ussror and STANLEY Corrsin 


SUMMARY 


An experimental and analytical study has been made of some 
Jeatures of the turbulent heat diffusion behind a line heated wire 
stretched perpendicular to a flowing isotropic turbulence. The 
mean temperature distributions have been measured with sys- 
tematic variations in wind speed, size of turbulence-producing 
grid, and downstream location of heat source. The nature of 
the temperature fluctuation field has been studied. 

A comparison of Lagrangian and Eulerian analyses for dif- 
fusion in a nondecaying turbulence yields an expression for 
turbulent-heat-transfer coefficient in terms of turbulence velocity 
and a Lagrangian “secale.” . 

The ratio of Eulerian to Lagrangian microscale has been de- 
termined theoretically by generalization of a result of Heisenberg 
and, with arbitrary constants taken from independent sources, 
shows rough agreement with experimental results. 

A convenient form has been deduced for the criterion of inter- 
changeability of instantaneous space and time derivatives in a 
flowing turbulence. 

INTRODUCTION 


One of the most striking aspects of turbulent motion in 
fluids is its dispersive property. This “convective diffusion,” 
illustrated by the general statistical tendency of (noncon- 
tiguous) fluid elements to get farther apart with increasing 
time, was probably first observed long before the era of 
, analytical fluid mechanics. An analytical start on this 
problem was not made, however, until the now-classic work 
by Taylor in 1921 on diffusion by continuous movements 
(reference 1). Not only did this paper lay a groundwork for 
the study of turbulent diffusion but it also represented a 
forward step in the ideas essential to development of a gen- 
eral statistical theory of turbulence, a field which had scarcely 
progressed since Reynolds’ original formulation of the equa- 
tions of motion for a flow in which mean and fluctuating 
parts could be distinguished. 

The diffusive action of a turbulent flow may manifest itself 
in various ways, depending upon the initial and/or boundary 
conditions and upon the interests of the observer. The fol- 
lowing possible measures of the diffusive powers are neither 
exhaustive nor mutually independent: 

(1) The average rate of dispersion of particles from a fixed 

source 

(2) The average rate of increase of spacing between differ- 

ent particles 


! Supersedes NACA TN 2710, “Diffusion of Heat From a Line Source in Isotropic Turbu- 
lenea”’ by Mahinder 8. Uberoi and Stanley Corrsin, 1952. 


(3) The average rate of transport of particle concentration 

under a given mean concentration gradient 

(4) The average rate of increase of the length of a fluid line 

(5) The average rate of increasé of the area of a fluid sur- 

face 
The word “particles” means simply indelibly tagged fluid 
elements, much smaller than the smallest length associated 
with the turbulence. 

The present report is concerned primarily with measure 
(1). The measurements have all been made in the thermal 
wake of a long thin heated wire mounted perpendicular to 
an isotropic turbulent air flow and producing no turbulent 
wake. Here the tagging is thermal, and the degree of in- 
delibility (negligibility of molecular diffusion) is one of the 
matters to be investigated. l 

The diffusive property (for scalars) of a turbulent flow is 
apparently a secondary characteristic at least in the sense 
thatit need not explicitly enter the dynamical problem. The 
diffusion may be regarded as a kinematic phenomenon, to be 
deduced from the dynamical solution. to the problem if and 
when the latter is obtained. Thus the objective of research 
on turbulent diffusion may be to seek a connection between 
the diffusive and the dynamical statistical variables, even 
before the complete dynamical theory is available. 

Measure (3) is usually termed the ‘turbulent transport” 


or “transfer” problem. Although of extensive practical im- 


portance, it has not yet been subjected to genuine theoretical 
study. 

Most of thesemiempirical “theories” of turbulent transport, 
for both scalar and vector properties, employ an Eulerian 
formulation of the basic equations, and up to now they have 
been unable to relate the turbulent transport correlation to 
other statistical . functions describing the flow. Taylor 
(reference 1) showed that in the simple case of a homogeneous 
field of isotropic turbulence, and even in a decaying isotropic 
turbulence (reference 2), a Lagrangian formulation of the 
transport (i. e., diffusion; the terms will be used interchange- 
ably) problem leads to some important results.? 

Up to the present time little theoretical or experimental 
work has been done to find relations, if any, between the 
Lagrangian statistical measures of a turbulent field with its 
Eulerian statistical measures. Since turbulence dynamics 
seems best handled in the latter terms and turbulent diffusion 


-in the former, it is evident that such a connection is impor- 


3In his tensorial generalization of Taylor’s work on case (1), Batchelor (reference 3) has 
chosen to call this an “Eulerian” analysts, describing only case (2) as “Lagrangian.” In 
keeping with previously accepted nomenclature, both cases are Lagrangian, (1) involvinga 
single particle and (2) involvinga pair. In fact, case (2) might be termed a mixed (Eulerian 
and Lagrangian) problem. í 


815 


ee a aa e —_—— 


816 


tant. Hence, one of the purposes of the present experiments 
has been to compare the magnitudes of some of these quanti- 
ties under variations in the turbulent field. For example, 
the postulates of Taylor and Heisenberg on a relation between 
Lagrangian and Eulerian microscales can be examined and, 
in corrected form, compared with experiment. 

The turbulent diffusion from a fixed line source can be set 
up analytically as an ordinary (Eulerian) “heat-transfer’’ 
problem, permitting a start to be made in relating measures 
(1) and (3) of the diffusive power of a turbulent flow, under 
certain simplifying assumptions. 

Measures (4) and (especially) (5) may well be classed as 
characteristic of the ‘‘turbulent-mixing”’ problem rather than 
diffusion in the common connotation. 

Experimental work on diffusion from a fixed local source 
in a turbulent flow has been meager. In isotropic turbulence, 
there have been the measurements of Schubauer (reference 
4), Simmons (reported by Taylor in reference 2), Dupuis 
(reported by Kampé de Feriet in reference 5), Frenkel (refer- 
ence 6), and Collis (reference 7). Of these, only the data of 
Simmons and Collis are extensive enough to permit confident 
computation of the Lagrangian correlation function.: In 
turbulent shear flow, Skramstad and Schubauer (reference 8), 
Dryden (reference 9), and the present authors (reference 10) 
have measured distributions close to a source; Kalinske and 
Pien (reference 11) and Van Driest (reference 12) have made 
measurements somewhat farther downstream. 

None of these studies was repeated with a systematic 
variation of the properties of the turbulence. In spite of 
the poor precision inherent in this type of measurement, it 
was hoped that such-an approach would at least show up 
some general trends in the relations between Eulerian and 
Lagrangian variables. 

This investigation has been conducted at the Aeronautics 
Department of the Johns Hopkins University under the 
sponsorship and with the financial assistance of the National 
Advisory Committee for Aeronautics. The authors would 
like to acknowledge the assistance of Messrs. Alan Kistler, 
George Stierhoff, and Allen Gates and Miss Patricia O’Brien, 
as well as-the helpful criticism of Dr. Francis H. Clauser and 
Dr. C. C. Lin. 


SYMBOLS 
c root-mean-square molecular velocity 
Cy ‘specific heat per unit volume at constant pres- 
sure 
f(r) Eulerian velocity correlation coefficient (nota- 
tion of Von Kármán and Howarth) 
He2pc,U f Ody 
0 
a 
pe, |. f (EdE 
3 width of rectangular heat pulse 
k thermal conductivity 
ker turbulent-heat-transfer coefficient 


L Eulerian scale (z = i f(r) dr) 
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Ly aerameanatale (t= Í “B(a)dr) 

L Lagrangian scale for nondecaying and decay- 
ing turbulence (1, = f Bda) 

l mixing length 

M grid mesh size 


Mim | Rian 


Ps) ‘probability density of fluctuating temper ature 
p static pressure 
Rz turbulence Reynolds number based on Euler- 
ian scale (R: == ee 
2 _ v@ot—7) 
tRi-r =RK,(, T)= v (Ov (t— T) 
R,(r) Lagrangian correlation coefficient for nonde- 
caying turbulence 
R(n) Lagrangian correlation as a function of n for 
nondecaying and decaying turbulence 
Ry =L 
6U 
Ry turbulence Reynolds number based on Euler- 


o'r 
lan microscale (m= > 


scalar distance between two points 


-GET 


8 average on-center spacing of pulses 
1/2 
r=|($ ae ir (a) | 
2) 

t 

U mean velocity in adron 

u instantaneous velocity fluctuation in z-direc- 
tion. 

au’ N ' 

Vv mean velocity in y-direction 

v instantaneous velocity fluctuation in y- 
direction 

v =D 

w instantaneous velocity fluctuation in zZz- 
direction 

X distance traveled in z-direction by a fluid 
particle i 

£ distance downstream from grid 

Lo location of heating wire 

Ar=2%—Zo 

Y distance traveled in y-direction by a fluid 
particle 

Y’ root-mean-square displacement of a fluid 
particle in y-direction (Y =Y") 

Va root-mean-square displacement of a molecule 

yY distance in direction of measured diffusion 

Z. distance in direction of heating wire 

v 
Oe +U Ox 
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(lef) = f v&d 


Q instantaneous temperature (measured above 
ambient room temperature) 


a] mean temperature 
0, maximum mean temperature at a cross sec- 
tion, a function of Ar 
B instantaneous temperature difference 
Vays 
Do temperature difference of rectangular heat 
pulse 
K dimensionless empirical constant 
A mean free path of a molecule 
A Eulerian microscale of turbulence 
1 a 
C=) ) 
Àr Lagrangian microscalé of turbulence for non- 
, 2 1/ 
decaying turbulence =(=) ) 
An Lagrangian microscale of turbulence for 
nondecaying and decaying turbulence 
22 1/ 
| C=( 0) ) 
y kinematic viscosity 
f= ria) ; \ 
p density 
time difference 
(7) mean value or ensemble average 
Subscripts: 
mar maximum 
min minimum 


EQUIPMENT AND PROCEDURE 
AERODYNAMIC EQUIPMENT 


The wind tunnel (fig. 1) is an open-return NPL type tunnel 
with a 2- by 2-foot working section and a free-stream turbu- 
lence level of v’ /U=0.06 percent and u’/U=0.05 percent at 
a mean velocity of 26 feet per second. The turbulence- 
producing grids were as follows: 


Dosignation 


l-in, grid.... | Biplane, bie ie 
in. -- | Woven, motal.. 
-in, grid... Woven, metal.. 


Square to ure transition sections-”” 


Over-all length = 425 ft 


““Cheesecloth screen 
All other screens ~ #24 wire mesh 
Figure 1.—Sketch of open-return wind tunnel. 
3216095 —05——52 i 
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They were mounted in turn at the upstream end of the work- 
ing section. 

The heat source was an 0.008-inch-diameter platinum wire 
stretched vertically across the tunnel at various distances 
from the grid. It was heated by direct current to tempera- 
tures between 500° and 700° C, with the latter figure only 
at the highest operating velocity of 38.0 feet per second. 
The wire Reynolds numbers at this condition and at the two 
other velocities were as follows: 


Reynolds number 
based on— 


A preliminary investigation was made without grids to 
insure that these operating conditions did not generate a 
vortex street downstream of the heated wire. 

With the grids in place, the mean momentum.wake became 
practically undetectable with total-head tube and manometer 
at distances greater than 1 or 2 inches downstream of the 
heating wire. 

MEASURING EQUIPMENT AND PROCEDURES 


The mean-temperature distributions were measured with a 
Chromel-Alumel thermocouple and a Leeds and Northrup 
type K-2 potentiometer. The cold junction was kept out- 
side of the wind tunnel. 

The shadowgraph technique was used to photograph the 
laminar thermal wake close to the source with no grid in the 
wind tunnel. This information was applied to the problem 
of “correcting” the thermal wake in turbulent flow for the 
effects of molecular diffusion and of finite source size. Figure 
2 (a) is a shadowgtaph of the wire wake with no grid in the 
tunnel; figure 2 (b) is a typical time exposure with grid- 
produced turbulence. A resistance-thermometer traverse of 
the laminar wake in a flow of very small turbulence showed 
that the temperature profile had already become very nearly 
Gaussian at a distance of 1 inch (125 wire diameters) down- 
stream. ‘The white lines on the sides of the dark wake 
shadow in figure 2 (a) correspond to the minimums in the 
second derivative of the density profile. For small tempera- 
ture differences these coincide with the maximums in the 
second derivative of the temperature profile. Although the 
temperature differences are not small in the immediate 
vicinity of the wire, this condition is reasonably well satisfied 
at relatively small values of Az as evidenced by the parabolic 
spread of this laminar wake. 

The standard deviation of the wake could then be com- 
puted from the spacing of these two bright lines, with the 
assumption of a Gaussian distribution. Since the closest 
points of traverse in the turbulent cases were \ inch (63 wire 
diameters) from the heat source, this was probably a reason- 
able assumption. 

As pointed out by Taylor Gerent 2), the molecular and 
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(a) No grid in tunnel. 
(b) Grid-produced turbulence. 
FIGURE 2.—Shadowgraph time exposure of wire wake. 


the turbulent diffusive phenomena are statistically inde- 
pendent, so that the squares of the standard deviations due 
to these two effects are additive. Hence the wake spread 
due to turbulence alone, from a true line source, was obtained 
by subtracting the square of the standard deviation of the 
laminar wake (computed from the shadowgraph) from the 
square of the standard deviation of the total wake (from 
thermocouple traverses with grid turbulence present) at all 
stations. This difference was the square of the standard 
deviation of the desired phenomena. Ail wake-spread data 
presented in the next section have been corrected in this 
fashion. 

Parenthetically, it should be remarked that the laminar 
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wake in figure 2 (a) spreads parabolically within the limits 
of precision, from at least 1 inch on, so that the effects of 
density differences on the flow phenomena must have been 
negligible for this investigation. 

The transverse turbulence levels v'/ŭ behind the grids 
were obtained from the initial rate of spread of the mean 
thermal wake (method due to Schubauer (reference 4)) after 
the effects of molecular spread and finite source had been 
removed. The resulting levels were somewhat higher than 
those obtained with a hot-wire anemometer but were used 
because of their consistency with the rest of the measured 
diffusion curve. 

Free-stream velocity fluctuations (without grids) and the 
wake temperature fluctuations (with grids) were measured 
with the hot-wire anemometry equipment described in refer- 
ence 10. The wires were 0.00025-inch platinum etched 
from Wollaston; the compensated response of the system 
was flat within +2 percent over a frequency range from 3 
to 12,000 cycles per second. 

Oscillograms of the temperature fluctuations were recorded 
by photographing a blue oscilloscope tube with fast 35- 
millimeter film in a General Radio type 651-AE camera. 

Probability densities of the temperature fluctuations at 
fixed points in the mean thermal wake were determined 
from photodensitometer traverses of time-exposure photo- 
graphs of a short-persistence (0.001 sec) blue oscilloscope 
tube with the temperature fluctuations on one pair of plates 
and a 30,000-cycle-per-second sweep on the opposite plates. 
The technique is essentially that used by Simmons and 
Salter (reference 13). 


EXPERIMENTAL RESULTS 
MEAN THERMAL WAKE. 


Complete mean-temperature wakes behind a line source 
of heat were measured for 10 différent conditions. Arranged 
to indicate the systematic variation of one parameter at a 
time, these conditions were as follows: 


sfo an ge 
Lin. grid; $7=43.4; Ü, fps l 


` Üm 25.6 {ps; <4 =43.4; grid, in 


* Um25.6 fps; == 86.1; grid, in 


U=25.6 fps; %-in. grid; TyM- 


Um 25.6 (ps; 4-in. grid; z./MM_....-. 43. 
72, 
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8, arbitrary scale 


ems 


Ə, arbitrary scale 


—, 


(a) a= 15, “*=43.4, M=1 inch, and U=26.6 feet per second. 


AY a, Zo i: F 
(b) VWs ii 7786.1, M Z inch, and U=25.6 feet per second. 
Frayre 3.—Experimental scatter; temperature behind line source of 
heat. 


Here U is the mean velocity, M is the grid mesh size, and 
a, is the heat-source location measured from, the grid. Since 
some individual cases enter as elements in two sequences, 
the total number of elements is more than 10. 

Two of the many mean-temperature traverses in the 
y-direction (perpendicular to mean flow and to source line) 
are shown in figure 3 to give an idea of the amount of exper- 
imental scatter. The upper traverse was the worst of the 
lot, even showing an apparent skewness which was not 
borne out by the investigation as a whole. Theower trav- 
erse is more nearly typical of the measured temperature dis- 
tributions from which the standard deviations of the mean 
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FIGURE 4.—Spread of heat from a line source. a= 48-4, M=1 inch. 
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FIGURES 5.—Spread of heat from a line source. apa 43.4, U=25.6 feet 


per second. 


thermal wake were computed. By comparison with the 
reference curve, it is seen to be essentially Gaussian. This 
was the case for temperature profiles at all stations. Since 
the virtually Gaussian character of such a wake has already 
been established by several of the-earlier publications, there 
seemed to be no point in reproducing here all of the large 
number of traverses measured. 

The mean-thermal-wake spread for the 10 different con- 
figurations studied is given in figures 4 to 7 as plots of cor- 
rected standard deviation Y’ against distance from the heat 
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FIGURE 6.—Spread of heat from a line source. y= 86.1, U=25.6 feet 
per second. 
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Figurp 7.—Spread of heat from a line source. 
feet per second. 


source Az. Each point in these figures corresponds to a 
complete transverse temperature traverse. 

In order to have values of transverse turbulence level »’/U 
consistent with the thermal-wake behavior, these values 
were determined from the initial angle of spread of the cor- 
rected wake standard deviation (references 2 and 4) instead 
of from direct hot-wire anemometer measurements. The 
results are plotted in figure 8. Since these represent an 
insufficient number of points per grid to permit the drawing 
of reliable curves, some simplifying assumptions were made 
based upon the results of several more-detailed turbulence- 
decay investigations. (See references 14 to 20. The perti- 
nent results of references 14 to 16 are summarized in refer- 
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FicurE 8.—Decay of turbulence as determined by thermal wake of a 
line source. 


ence 17.) In the light of these papers, it was assumed that 
the decay curves had a common apparent origin and that 
this was obtainable by drawing the best straight line for all 


- the available points, independent of the differing wind speeds 


and mesh sizes. In computing any individual Lagrangian 
correlation function from the corresponding wake history, 
the turbulence decay rate was assumed to be given by the 
line drawn through this common origin and the specific tur- 


` bulence value giving the measured initial spread angle for 


this wake. This is, of course, a very rough procedure, but 
the experimental scatter in this whole method of determin- 
ing Lagrangian correlation functions is so great that a more 
extensive study of decay (including the resolution. of. incon- 
sistencies between wake method and hot-wire method) 
seemed unwarranted at this time. 


TEMPERATURE FLUCTUATIONS 


Distributions of temperature-fluctuation level #’/6 in the 
thermal wake have been measured by using the hot-wire 
anemometer as a resistance thermometer, that is, at a 
current low enough to render the sensitivity to velocity 
fluctuations negligible compared with the sensitivity to tem- 
perature fluctuations (reference 21). A representative dis- 
tribution of #’/6 in the z-direction is given in figure 9. Typ- 
ical transverse distributions are given in figures 10 (a) and 
10 (b). Itis clear that the temperature-fluctuation inten- 
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O 10 20 30 40 50 60 TO 


Ax, in. 
Fiaur5 9,—Temperature fluctuations at 7,43 inches, M=1 inch, and 
U=25.6 feet per second along wake axis. 


80 90 


bat 
Ə 


T a 


(a) Az=10 inches. 
IiaurE 10,—Temperature fluctuations behind a line heat source at 


3743, M=1 inch, and U=25.6 feet per second. 


sity changes very little with increasing values of Ar. A 
rough explanation of the very high values of W/O (com- 
pared with the concomitant turbulence level, for example) 
in terms of the highly intermittent structure of the thermal 
wake has been given in reference 10 and will be discussed 
in more detail later in this report. This intermittency is 
shown very clearly in figure 11, a series of temperature 
oscillograms recorded at two different positions across the 
thermal wake for a fixed value of Ar and at two different 
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Figur 11.—Oscillogram records and normalized probability density 
of temperature fluctuations. Film speed, 1 foot per second. 


values of Az with y=0. The one-sided and pulse character 
of the instantaneous temperature at a fixed point in space is 
also demonstrated by its probability density. 


THEORETICAL CONSIDERATIONS 


HOMOGENEOUS STEADY TURBULENCE AT REST 


For a nondecaying incompressible turbulence with no 
mean motion, Taylor (reference 1) was followed in getting 
an expression for the mean time rate of diffusion in the y- 
direction (say) from a fixed source as measured by the 
second moment of the probability density of the diffusion, 


that is, the mean-square particle displacement Y(2): 


929 
da aF 
= 200 f ' vdt 


The bar denotes ensemble average. 

Taking v(t) inside the integral, interchanging the processes 
of integration and averaging, and introducing the Lagrangian 
(auto) correlation coefficient (following the fluid particle), 
where r=i—t, 


Rijs v(t) v(t— 7) 
f v? 
there results l 
ESS t PT 
Fi=27 f f R,()drdT (1) 
0/70 


This is Taylor’s form. Integration by parts yields a form 
like that in the work of Kampé de Feriet (reference 5): 


Yi= 20 f (¢— T)L(r)dr (2) 


In this Lagrangian analysis, v(t) is the velocity of a fluid 
particle in the y-direction at time t; »({—7) is the velocity of 
the same particle at time ‘—r. Corresponding expressions 
can be written for the rate of diffusion in any direction. 

Diffusion from an infinite line source, the case to be 
discussed here, is a two-dimensional problem in the mean, 
and, in addition to equation (2), 


X= Qu? f (t—rT)Eulr)dr (3) 


A tensorial generalization of these concepts has been given 
by Batchelor (reference 3) dealing with the behavior of 
AXiA,(t) where X, and X; are any two of the orthogonal 
displacements of the particle at time t. 

It should be noted that this analysis gives no information 
on the shape of the probability density of Y(é) or of X(). 
In fact there still exists no theory for these. However, 
experiments in flowing turbulence (the case to be considered 
next) show Gaussian density, within the PEA 
precision, for Y(t) at all values of è. 

For this stationary random process, Taylor (reference 1) 
introduced the.concept of the Lagrangian “scale,”’ 


io Í "Reji ac 


These have the dimensions of time and are characteristic 
constants of the system. 

In his later work on the (Eulerian) dynamics (reference 2), 
Taylor had occasion to introduce another measure of the 
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correlation function, which he called the ‘‘microscale.’’ 
Applying the same geometrical concept to the present 
function, the Lagrangian microscale 


2 
EO), 3 2 


is simply the r-intercept of the vertex-osculating parabola 
of the even function &,(r). The kinematic significance is 
clearly shown by a series expansion of v(t+r) in R,(r): 


Aro} oo+| (w) fet m 3) 
GE CaN! Hae (ai), fet | 


But = Constant, so that 


v(t)v(t-+7) 
v? 


Ro)=1—-5 l LEN pais (6) 


From equations (6) and (5), the Lagrangian microscale for 
v(t) is 


~ 


- 2 
Nima (7) 
a) 
or = 
d 2 2 
a =2 (Ta) 


TURBULENCE IN A FLOWING MEDIUM 


The dictates of both practical interest and experimental 
feasibility require analysis of the diffusion when there is a 
mean. velocity U relative to the source. Since the diffusion 
phenomenon is linear, the probability density (mean-concen- 
tration distribution of tagged particles in the wake) is simply 
proportional to the superimposed probability densities of a 
continuous line of sources moving with the mean velocity U 
with their time (and space) origin at e F fixed source. 


This is illustrated in figure 12(a) for = 7 ara The circles 


(corresponding to: isotropy) are the standard deviations of 
the dispersions that would occur from moving sources. The 
envelope of these circles gives a measure of the mean wake. 
It is obvious that in general the functional form of the mean- 
concentration distribution along a line Ar=Constant will not 
be the same as the functional form of the same’ quantity for 
the individual source at time ¢, that is, at position z=U¢. 
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Fraure 12.—Propagation of turbulence from source moving with mean 
velocity U. 


/ 
However, it is extremely likely that if a <i (fig. 12(b)) 
the mean concentration along a line Ar=Constant becomes 


very nearly of the same functional form as that for the source 


at a=Ut. For a zero-correlation Gaussian density, the 


equivalence is easily demonstrable.* 
/ 


The condition a <i will always occur at large enough 


values of Az when F<. This follows from the asymptotic 


3 The most gencral mathematical restrictions under which the superposition of a line of 
identical densities will yield a “cross-section” density of the same form have not been studied 
here. It & obvious that the condition of statistical independence ts sufficient. The fact that 
the density of each of the velocity components in isotropic turbulence has been found to be 
Gaussian within the experimental precision’ seems to show that the equivalence under 
discusston is at least a good approximation. 
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parabolic behavior, Y” (¢) ~t (reference 1), of the diffusive 


process. It will occur for all values of Az when F<. There- 


fore, in this case a simple approximate space-time transforma- 
tion în the mean is permissible, and the ¢-variation in Taylor’s 
theory of diffusion by continuous movements becomes a 
variation of z/UY. (It must be emphasized that the fore- 
going discussion does not apply directly to the possibility of 
applying a space-time transformation to the instantaneous 
turbulence variables. This latter question will be discussed 
later.) 


z 
In the present measurements F< and, therefore, the Ar- 


variation of diffusion gives an approximate measure of the 
Lagrangian correlation coefficient (in time). Equation (2) 
can be rewritten as: 


V(x mo g fe- - © 


DECAYING ISOTROPIC TURBULENCE 


When the turbulence is decaying in time (similar to space 
in the flowing turbulence) v? is no longer constant, and the 


analysis cannot be carried out as far as equation (1). The 
same approach stops with 
 — t 
GV =0'O [ot DR (9) 
t 0 


where the prime denotes root-mean-square value, and 


v(t)vg—r) 


Bet 9100 1) 


There is no a priori reason to believe that ,#,_, is a function 
of 7 alone, as in the nondecaying turbulent flow. 

At this point Taylor (reference 2) invokes the empirical 
fact that over a wide range of mean velocities (all of which 
give essentially the same distribution of »’/U in x behind a 
grid) the thermal wake behind a line heat source at fixed x 
appears to be unchanged in form, within the experimental 
error. This is consistent with dependence of the diffusive 
process upon a variable of type 


f 


t 
es Í, v't) dt 


=gh" ik ax, 


Therefore, Taylor has postulated the unique dependence of 
LR, on the variable n. With this postulate and the space- 
time transformation valid for small turbulence level, he 


(10) 


824 


arrives at 
Fe P2 S| R) dm (11) 
or ` 
oy °R d 12 
dat 3(71) 71 ( ) 
where 
d zu, 
7 U an 
=v'dt 


Equations (11) and (12) look like the equation for non- 
decaying turbulence. They also give 


Ying f “oa a (13) 


A physical significance of the length 7 is underscored by 
the limiting form of equation (13) as „—>0 and R,—>1. Then, 


Y’=n 


so that n is a measure of the lateral diffusion that would 
occur if the lateral velocity fluctuation following a particle 
v(t) remained perfectly correlated but decreased in magni- 
tude according to the decay rate of the turbulence level. 

In postulating ,2,.=R,(n), Taylor was apparently com- 
paring g diffusive processes in turbulence fields with 


identical = iii (x). Of more general interest is the comparison 


of diffusion in fields with differing turbulence-level distribu- 
tions. Although such a generalized application of his postu- 
late is doubtless not too well applicable, it is conceivable it 
might have. approximate success in the more general com- 
parison. For both convenience and lack of any obviously 
superior alternative, his suggestion is therefore applied in 
computing the results of the measurements reported here. 4 
With the y-postulate, the R, Lagrangian correlation func- 


tion can be obtained from measurements of Y? as a function 
of Az: 


Y? (14) 


t After this work was compkted, Dr. Batchelor suggested an alternative approximate 
approach: In order to construct a stationary random function out of the nonstationary v(t) 
one first normalizes the dependent variable with its root-mean-square value. (This has been 
nutomatically accomplished by use of the correlation coefficient.) If the decaying quantity 
is assumed to maintain complete similarity during decay, all characteristic times (e. g., 

time scale and microscale) vary in the same way with ¢ and the new independ- 
ent variable is constructed by dividing r by this ¢-variation. Unfortunately, this variation 
is unknown a priori, so it would be necessary to assume further that the Lagrangian scales are 
directly proportional to the Eulerian scales, about which there is previous experimental infor- 
mation. His work has now been published; see reference 22. 
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A scale and a microscale can also be defined for R;,: 


L= | Rina (15) 
2 __ 2 
Ae a R0)” (1 6) 


With nondecaying turbulence, no one-parameter true 
(time) Lagrangian correlation function exists, and the 
n-formulation is much more convenient. A further signifi- 
cance of this variable will appear in the comparison of 
Eulerian and Lagrangian treatments of diffusion from a line 
source in flowing turbulence. 


ACCELERATIONS IN DECAYING TURBULENCE 


A series expansion of v(t-+7) for decaying turbulence will 
show something about the initial behavior of the true (time) 
Lagrangian correlation function and will indicate an experi- 
mental method for examining a hypothesis of Taylor on the 


_ interchangeability of instantaneous time and space deriva- 


tives when the turbulence level is low (reference 23). 
Write the Lagrangian correlation coefficient 


_ »Ort—7) 
R(t, r)= v (Àv @— 7) 


a) +(e) at” 


into numerator and denominator, and restrict the analysis 
to small values of 7: 


MOREEN 


Divide numerator and denominator by v4, expand the square 
root in the numerator, and keep terms in 7°: 


| Benmi- (F) H(t), ba 


For negligible decay rate this reduces to equation (6). 


(17) 
Substitute 
vi—r)=v( 


£,(t,7) (18) 


(19) 


Equation (19) shows that a Lagrangian microscale defined 


7 Ry(t, 7) 
1 1— all, T 
rO A im | =e =en] 
oF, (t, T) 
-| N] (20) 
is expressible as 
1 1 do 
sora G- #) | a) 
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Introduction of Taylor’s y-postulate transforms equation 


(19) to 
| dv\? a 
F(a) er tale ya o A 202_, 


after the additional approximation that v’(n) =v (0) when 
n is very small. 
This gives a new expression for Lagrangian microscale ^: 


oa HE dvy\? 

dy (T) 
Equation (23) is in contradiction to. equation (17) of part 
IV of reference 2. In that work Taylor has apparently 
assumed that v’ is a nondecaying function of n. However it 
certainly is decaying, even in terms of this distorted coordi- 

nate, and the dv’/dy-term must be included. 

Since A, and v’(z) can be determined experimentally from 
the mean thermal wake behind a line heat source, equation 


(22) 


(23) 


PANS 
(23) permits determination of Ẹ which is simply related 


to the mean-square “Stokes” acceleration GO} J A =o (Z N. 


This quantity is of particular interest for the possibility of 
an instantaneous space-time transformation. at low turbu- 
lence levels. This was first proposed by Taylor (reference 23) 
and has since been used very widely, especially to get approx- 
imate values of partial derivatives with respect to x (the 
mean-flow direction) by measurements of time partial 
derivatives. 

The total (or Stokes) derivative of v(z, y, z, ġ in a turbulent 
flow with mean velocity U along the z-direction is 


a 


__ Ov on, dv 
= S2 LT) S2 Togt oy (24) 


Taylor’s hypothesis amounts to the statement that 


_ 


or 


with 
Ov? 
U (dz 
In detail, equation (25) is 


I Ov 1/2 
tdt Ox; 


Te. 
“bua F T 


P (3) 


«&I1 (25a) 


(24a) l 
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In the absence of information on the algebraic sign of the 
triple correlation term, it is sufficient to require the two 
conditions.* 


’ (26) 


can be determined from measurements of Y’ (x); 


But (2) 


JON? Lo 
or = 2a 
and upper bounds, in terms of measurable functions, can be 


where A is the Eulerian microscale (reference 2); 


set on uu; eee with the use of Schwarz inequalities. Thus, 
Ox; Ox, 


an experimental check of the requirements in equation (26) 
is to be made in the section entitled “Computation of 
Results.” 


RELATION BETWEEN EULERIAN AND LAGRANGIAN MICROSCALES 


Taylor (reference 2) inferred an approximate relation 
between à, and à by neglecting the effect of viscosity on 


FIE a 
pressure-gradient fluctuations and estimating [E S2) | es 


being approximately Bele (22) Te This led to a constant 


ratio A,/A for all turbulence. The rough nature of this 
analysis induced Heisenberg (reference 25) to conduct a 
more detailed study of the static-pressure fluctuations and ‘ 
to reestimate the d,/A ratio. However, he followed Taylor 
in ignoring the dv’/dy-term in the relation between i, and 


. 2 
(F (see equation (23)) and in neglecting viscous terms in 


the relation between ^, and (3 


Although these omissions are probably not serious except 
in the low Reynolds number range, it seems interesting, if 
only for the sake of completeness, to use a Heisenberg type 

———— . 
of approximation for (52) and to repeat his treatment with 


the omissions rectified. 


+ Lin has discussed the validity of Taylor’s hypothesis using a slightly different formulation 
in reference 24. He points out there that if A= È as, then (from tho Schwarz inequality 


(4) = È qa” TEJ where a, Is a set of numbers, 


$26 


From the complete Navier-Stokes equations in the wave- 
number space, Heisenberg deduced an approximate expres- 


* gion for (Vp)? in terms of mean quadruple products of the . 


“harmonics” of the velocity field. His principal simplifying 
assumptions were: 

(a) Different Fourier components of the velocity field are 
uncorrelated 

(b) The tyrbulent energy spectrum is given by the solution 
to his equilibrium-energy-transfer equation, above a lower 
cut-off wave number ke 

Following these, but using Chandrasekhar’s (reference 26) 
solution to the Heisenberg equation instead of the inter- 
polation formula used by Heisenberg, there results 


(op. _ 25. as 
(ioe 
where x is a dimensionless empirical constant. 

The numerical constant in equation (27) would perhaps 
have been given more accurately by the use of a “‘self- 
preserving” spectrum (calculated by Chandrasekhar from 
Heisenberg’s equation) instead of the stationary spectrum 
with low cut-off wave number. Time was not taken to make 
the requisite additional calculations because: (a) The value 


2, 0Y 
a (27) 


AmA 2 
of (5 depends principally upon the high-wave-number 


region of the velocity spectrum rather than the low-wave- 
number region, where the difference would be greatest, and 
(b) the experimental results and (especially) the value of 
x both have a considerable range of uncertainty. 

The mean square of the y-component of the Navier-Stokes 
equation will lead to a relation between ñ^ and A: 


oe oe Ov Ov 
ot ee 


therefore, 


where the correlation between pressure gradient and velocity 
Laplacian function is zero because of isotropy. 

Equations (23) and (27) give the first two terms in terms 
of the microscales, and the mean-square Laplacian function 
is expressible in terms of the fourth derivative of the Von 
Kérmén-Howarth f(r) correlation coefficient at r=0 (refer- 
ence 27): 


Ways Vs") 29) 


Consequently, equation (28) becomes 


(vt | (dv’\?_ 25.2 (vF 
eas x)= K 


+2 (oF) 


The second term in this equation can be replaced by the 


(2) =; (32 (22) HAr (28) 


x=0.45 has been used here. 
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turbulence-decay equation 


d(v’}* (v’)" 
whence 
do'e v’ 
dt at 
so that 
2 oF Sort 25» 2 ey _ 25. 2 224 == pelona fto) 
therefore, 
Ae T260 1 12535 aja 
E r Ro Re LEA J= (0) Ca (31) 


é 
vr 
where hy=— 


Batchelor and Townsend (reference 19) have deduced an 
expression for f*"(0) which is valid in the region of decay 
where both 1/(v’)? and X increase linearly with ¢ (correspond- 
ing to large values of Ay): $ 


MPO= +s BS (32) 


where S=— y (2) / | (22) D) n » the skewness factor. Their 


experimental results showed S=0.39 approximately constant 
for isotropic turbulence. Then the estimate for MM 


becomes 
(G)-¥ Paes mol 4114 ie (83) 


The value of x was first estimated by Heisenberg (reference 
25), from measurements of turbulence decay, as 0.85. 
This method may be regarded as emphasizing the (relatively 
low wave number) energy-bearing range of the spectrum. 
Lee (reference 28) worked out an estimate based upon 
skewness factor (x=0.13), which gives heavy weight to the 
high-wave-number range. Proudman (reference 29) has 
reestimated x by comparison with measured curves of the 
double and triple velocity correlations. The value .=0.45 
leads to reasonably good agreement for the moderately 
high-wave-number region, over a wide range of values of Ry. 

It may be remarked that the supposed constancy of x is 
merely a postulate of the Heisenberg dimensional formulation 
of the spectral transfer function. -In fact if is by no means 
obvious that this turbulent part of the transfer is quantita- 
tively independent of the amount of spectrally local dis- 
sipation to heat. In any case, Proudman’s estimate of 
Therefore, 


pas 


Oy 12.5 7 (34) 


In the limit of Ry — 0, sudan (34) does not apply since 
equation (32) does not apply. However, the appropriate 
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limiting relation can be obtained directly. In this limiting 
condition the pressure term in equation (28) is negligible 
compared with the viscous term (the former goes with 
1/R, and the latter with 1/R,?), and the Eulerian velocity 


correlation coefficient is (reference 27) f(r)=e7"/™’. This 
gives f" (0)=3/M and 
A\* 3.5 
($) =z (35) 


EULERIAN ANALYSIS OF HEAT DIFFUSION FROM A LINE SOURCE 


The two-dimensional turbulent-heat-transfer equation is 


=00_k + (0 a} 
Oy pc,\oz? dy? 


Z o) 


where 6 is mean temperature, 3 is temperature fluctuation 
about the mean, k is thermal conductivity, and c, is specific 
heat at constant pressure. 

For the thermal wake behind a line source in isotropic 
turbulence with constant mean velocity V=0. With 
restriction to low turbulence level, a “boundary-layer” type 
of approximation can be applied to the mean wake, so’ that 


5z (ou) — (36) 


and 


so that equation (36) takes the approximate form 


00k 0% lige 
ee ply oy? dy pe 


(37) 

It must be emphasized that for this particular initial con- 
dition on the temperature (effectively a “point source”), 
the restriction to small turbulence level »’/U<i does not 
imply that 9’/6 is small. In fact, for this problem ¥#’/0 is 
often greater than unity, especially at the “edge” of the mean 
wake, as has been discussed in reference 10 and will be 
brought out again later in the present report. 

When the molecular transport can be neglected relative to 
turbulent transport 


7,00 ð Fo) 
Uo y Üy 


(38) 
an equation given in reference 10; a slightly more general 
treatment follows: 

With a constant rate of heat generation (similar to steady 
state in the average), the application of a Von Kármán 
integral-relation treatment to equation (38) yields an integral 
condition: 


2pc Tf © dy=Constant=H (39) 


where H is the average time rate at which heat crosses all 


planes perpendicular to U per unit length of heat source - 
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(z-direction). Of course, du has been neglected relative to 
OU in equation (39). 

Equation (38) has two unknowns, and the first objective 
is to express ® as a function of the (more easily measurable) 
Oley). After integration with respect to y, 


— >10 
vo=—U R ap Yt) 


But, by symmetry, 3v=0 for y=0. Hence F(z)=0 and 
(40) 


This relation is sufficient for the computation of #v(z,) 
from the measured Olzy) but the empirical fact of simple 
geometrical similarity in ©(z,y) suggests exploitation of the. 
consequent simplification. 


Assume 
Q(z, y)=8,(z) SE) (41) 
where §=y/Y’(x). This transforms equation (39) to 
= H* 
TONA (42) 


where H*= =Constant. It transforms equa- 


2ps | © dë 
tion (40) to 


Fo=(— U Tf sO d) Ti) G Be (T eg EJ E BD T 


(43) 


With equation (42), 6,(x) can be eliminated from equation 
(43), and after integration of the dimensionless integrals this 
leads, to the final form for the turbulent-heat-transfer 
correlation, 


1 dY’ 


vlz, Y=H* tf) == VG ar (44) 


where 


en Se 
i= Y"(xr) 


The same sort of analysis can be made on equation (37) 
which includes the molecular conduction, but the rather large 
experimental scatter in the present measurements seems to 
make such a refinement inappropriate. 


An “exchange” coefficient or “diffusion” coefficient for 


turbulent heat transfer kr is simply expressible in terms of 


Q(z, y). A conventional procedure for semiempirical analyses 
is to write for the turbulent transport an expression just like 
that for the molecular transport: 


— dO 
Sote E (45) 


which serves as the definition of kr. 
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For the simple case of equation (44), it turns out that 
— dy’ 
bpt a Y' (e) = (46) 
dé 


This has the particularly interesting property that in a 
nondecaying turbulence at very large values of z, where the 
mean thermal wake spreads parabolically (reference 1), kr 
becomes independent of explicit dependence on z. 

A more startling simplification follows for a particular 
. mean temperature distribution across the wake: All depend- 
ence of kr on y disappears if 


af 5 
de ef 
that is, if 
f= (47) 


But this is the Gaussian function, which is found empirically- 


to fit all the measurements within the experimental scatter. 
Hence one arrives at the empirical result that in both non- 
decaying and decaying turbulence kr is independent of y 


in the thermal wake behind a, line source of heat. From 
aquntion G7), Ser and 

f EE D 

kr=ptU Y (2) ae. (48) 


It can be seen that in the nondecaying case at very large 
values of z, kr is constant and independent of both y and z. 
RELATION BETWEEN SOME LAGRANGIAN AND EULERIAN 
PARAMETERS IN TRANSPORT 

There has apparently been little effort to relate the 
Eulerian and Lagrangian formulations of turbulent diffusion 
up to the present time. Exceedingly simple boundary con- 
ditions permit some connection to be made in restricted 
ranges of the present problem. 

For nondecaying or (with much less accuracy) decaying 
turbulence, equation (13) applies: 


=f (n — n) Palm) dn: 


or, equation (12), 


dy? , dY" 
d ~ dy 
=2 |" Radda 


z 
For small values of 7, R= 1 =a whence 
y 


Y?= 1? (1—5 5) (49) 
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dy? _ 1 7° 
PA i 


- Substituted into equation (44) these ” 


m=O% 4G- P, 2 (51) 
PEE (46) becomes 
tr=p00'y (1-3) (52) 
Janea f E T E et (1 p 1) In fact 
U Jo Y'(x) 7 12 2,? 


~ p’ 
for small values of 7, 7 =F aA 


/ 
For low-level nondecaying turbulence, =F x, and equa- 


tion (51) becomes 


One Y _1 vt af U 
a r E (53) 
while equation. (52) becomes 
O fp 
kom e0'y arias (1-a i) oa 


As 7-0 both equations (51) and (53) reduce to 


Fo= Hs) = 4, (55) 
while equations (52) and (54) reduce to 
kr= PCyV ay (5 6) 


At the other extreme, when 7 is very large, 
V?=2nL,—2M, (57) 


where Mi= f 7f,(n) dn =Constant and 


dY? 
a = 2L, (58) 
In this case, 
yL, a 
39 =H* fE) [2(n L, — M)? (59) 
and 


yL, 


O JOO y 
kr= pepe" ride BGL MI” (60) 


As indicated previously, the 7-variation can be expressed 
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in terms of z or t. 
If 7 is allowed to become large enough to make yL;> M, 


v’ 
_ YT 
dv =A" SO) sr RTP (61) 
1/2 
k= pepo! ey z) (62) 


The above formulas take on particularly simple forms if 
the empirical result of a Gaussian f(£) is utilized: 


J= 


Then the general expressions for dv and kr (equations (44) 
and (46)) become 


H* P 
volz, y)= EFO ay eal y _ (63) 


and 
na dY = 
kr(z,y)=00, Y0) —U (64) 
the latter having been deduced in the previous section. The 
particular forms for small values of n would follow from sub- 
stitution of equations (49) and (50) into these two. 


However, the most interesting form occurs for very large 
values of 7, There Y'(x) is given by equation (57) and 


H*yL,= ae 
volz, ) eens (65) 
Y= BG@L, My? OP GLM 


kr(z, y)=pc,v’ L,=Constant (66) 
For still larger values of 7, such that M, &nLy, 


aà = 


y 
n= 


and 
krz, y)= pcv L, =Constant 


The constancy of kr for large values of n (large values of t 
or x) is to be expected; a treatment of molecular diffusion by 
this method must certainly yield a constant coefficient for 
times much larger than the mean free time (of flight) of the 
molecules—that is, for all “macroscopic times.” Put an- 
other way, the simple parabolic behavior of Y’ for large 
values of 7 is a sure indication that © obeys the simple 
classical diffusion equation with constant coefficient, when 
viewed extremely ‘‘coarsely.”’ 

Perhaps the chief interest of equation (66) is its identifica- 
tion of L, as a significant Lagrangian length for diffusion at 
a large distance from the source. It enters the expression 


and Rl) . 


a * 
2/2 PE exp ( — itr) (67) 


for turbulent-diffusion coefficient in much the same way as 
mean free path enters the expressions for the molecular- 
diffusion coefficients. Furthermore, its role appears to be 
much like that attributed to Prandtl’s ‘‘mixing length,” 
which was brought into the turbulent-transport problem in 
a more or less intuitive fashion. 

Of course, the possible crude nature of Taylor’s original 
7-postulate may render the significance of L, more qualitative 
than quantitative in the case of decaying turbulence. 


COMPUTATION OF RESULTS 


Although Taylor’s assumption of the unique dependence 
of R- upon 7 is not likely to be accurate for collapsing 
together cases with widely differing turbulence decay rates, 
it does provide a relatively simple relation between Y’(zx) 
Therefore all of the mean-thermal-wake data 
were reduced on the 7-basis. i 

In principle tbe complete R,(n) curve can be obtained from 


Y?(x) by double differentiation (equation (12)): 


=! Ca 
or , ; 
1/UY & = 
Baal) ae 


However, simple double differentiation of the squares of 
a curve as uncertain as Y’(x) seems almost hopelessly in- 
determinate—although Taylor (reference 2) and Collis 
(reference 7) have apparently followed this procedure. A 
somewhat more circumspect technique has been tried here: 
The values of ^, and L, were determined first, through certain 
limit relations (to be described). Then the R, curve was 
determined by double differentiation, subject to the restric- 
tions of agreement with the previously determined scales. 
Thanks to rather poor determinancy of values of A, and L, 
this method is not so much of an improvement as it might 
first appear. 

LAGRANGIAN MICROSCALE M 

If equation (13) is restricted’ to very small values of 7 the 

parabolic approximation for R, can be introduced: 


reo (n— 11) Ce dm (69) 


therefore 


whence 


EEr (70) 


_ The computational procedure was to plot Y2/7? against 7? 
and to estimate the slope of the faired curve at n=0 where 
the curve must pass through unity. The abscissa intercept 
of ‘the 0-tangent is 6A}. The actual points, faired curves, 


7, in. 
(a) »>,=0.36. 
(b) L,=0.43. 
FiaurE 13.—Determination of M and ly = ar 43-4, M=1 inch, and 


U=8.5 feet per ee 


and tangents for all cases are presented in figures 13 to 22. 
Clearly the precision is poor. 


LAGRANGIAN SCALE L, 
Consider equation (12) in the limit as yo. Itimmedi- 
ately gives 
:=5 Tiii (71) 


2 y= 


and the graphical procedure based on this is also presented in 
figures 13 to 22. Some of the asymptotic slopes drawn are 
not the best representation of the experimental points. This 
is due to the auxiliary (assumed) restriction that R, cannot 
increase with increasing values of 7 as long as F, has not 
previously dropped below zero. The graphical precision 
attainable is perhaps a little hetter here than that for A}, but 


n, in. 


(a) M=0.34. 
(b) Ly=0.33. 


`- FiaurE 14.—Determination of A, and Lp; 7343.4, M=1 inch, and 


U=25.6 feet per second. 


the square root necessary to get A, means that A, is determined 
about as well as is L,. 


LAGRANGIAN CORRELATION FUNCTION R,(») 


With A, and L, determined, the initial (small Az) and final 
(large Az) behavior of the curve Y’(Az) is prescribed. ‘These 
parts of the curve were drawn on a graph witb the experi- 
mental points. Then the fairing in of a reasonable cen- 
tral portion to this mean Y’(Ar) curve was a relatively 
simple matter. The R, curve was then obtained by double 
differentiation. 

The curves drawn for Y’(Az) in figures 23 to 32 were de- 
termined in the fashion described above, as were the curves 
for R, in the same figures. 


EULERIAN MICROSCALE à 


In view of the approximate nature of the determination of 
^y, DO new direct measurements were made of the Eulerian , 
microscale A. Instead, \ was computed with tho energy 


DIFFUSION OF HEAT FROM A LINE SOURCE IN ISOTROPIC TURBULENCE 


7, in. 

(a) \y=0.27. 
(b) L,y=0.33. 

Fiaure 15.—Determination of Mẹ and Ly: a= 43.4, M=1 inch, and 


U=88 feet per second. 


equation for isotropic turbulence, from the measurements of 
turbulence decay: 


2 
Y= 10r (72) 
dv?/dt 
or, with the space-time transformation, 
ye 
N= — 10» —— (73) 
== dv? 
TEC 
dz 


EULERIAN SCALE L 


Earlier investigations have shown that the Eulerian scale 
in a grid-produced turbulence is closely a linear function of 
the mesh size of the grid producing the turbulence (for a 
given value of x and grid geometry) and is not significantly 
dependent upon the mean velocity (or grid Reynolds num- 
ber, provided it is sufficient to cause turbulence): Therefore 


the values of Z have been deduced from earlier measurements | 


S31 


7, in. 


(a) \,=0.25. 
(b) £,=0.34. 


Ficurs 16.—Determination of A, and Ly; So 13.4, M => inch, and 
U=25.6 feet per second. 


at the California Institute of Technology (reference 18) on 
grids of essentially the same geometry. _ 

Table I summarizes the results for Lagrangian and Eulerian 
scales and microscales. The results have been grouped to 
show the effect of systematic variation of one parameter at a 
time. Some of the results are presented in figures 33, 34, 
and 35. 


INSTANTANEOUS SPACE-TIME TRANSFORMATION 
The permissibility of an instantaneous space-time trans- 
formation in flowing turbulence, 


SF (74) 


can be estimated in accordance with’ equations (26). For 
equation (74) to be valid, the sufficient requirements are those 
given in equations (26) that: 
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A 


7, in. 
(a) 4=0.12. 
(b) £,=0.21. 
FIGURE 17.—Determination of >, and Ly; a= 43.4, M=} inch, and 
U=25.6 feet per second. 


avn 1/2 
i) 


ov\? 


With the aid of equation (23), the turbulence decay equa- 
Ap? ‘ne 
tion, and the Taylor relation ($2) = 2 +? the first of these 


(a) Ay=0.22. 
(b) L,=0.23. 


FicuRE 18.—Determination of M and Ly; 286.1, M=1 inch, and 
U=25.6 feet per second. 


conditions can. be written in the form 


O ARRO] o 


From the Schwarz inequality, essentially the necessity that 
the magnitude of any correlation coefficient be less than or 
equal to unity, 


Ov Ov ‘a 
“NITE, dz Uit =) (57) (76) 


where the prime in this expression denotes root-moan-squaro 
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min. 
(a) M =0.24. ; 
(b) D,=0.26. 
Ficurb 19.—Determination of , and La; =2=86.1, M=} inch, and 


U=25.6 feet per second. 


value. 
written 


For isotropic turbulence, equation (76) can be 


uny 92 98 Bo $25 (94473) SP (77) 


Thus, the second condition in equation (26) will be satis- 
fied if 


p’ 
2.7 = %1 78 
7 | (78) 


Both T and v'/U for the flows studied are presented in 
table I. It is clear that for these flows instantaneous gv and t 
partial derivatives may be taken proportional with reasonable 
confidence. 


833 


Po 
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(a) \,=0.11. 
(b) Ly=0.19. 
Fravan 20—Determination of M and L; 32=86.1, M=} inch, and 
U=25.6 feet per second. 


20o—CORRELATION 


The Eulerian measure of transverse turbulent heat trans- 
port is computed from the mean temperature distribution. 
The dimensionless form, #v/0,U, is given for two typical cross 
sections in figures 36 and 37. 


The measurements of v’/U and of 3/6, permit calcula- 


tion of the correlation coefficient Re=sen and this is also 


given in figures 36 and 37. 

For the convenient and reasonably accurate assumption 
of Gaussian mean temperature distribution, the correspond- 
ing turbulent-heat-transfer coefficient kr follows from equa- 
tion (64). It was found to*be independent of y, and typical 


O .05 10 45 ~ .20 30 30 ato) 40 
7, in. 
(a) ^= 0.22. 
(b) L,=0.22. 
Fieur® 21— Determination of M and Lr; Şe=172.3, M =5 inch, and 
U=25.6 feet per second. 


curves of kr/k are given in figure 38. The data for kr at 
three different speeds behind the l-inch grid are roughly 
collapsed together through division of kr by pcpv L, as 
suggested by equation (66a), an asymptotic result- for 
nondecaying turbulence (fig. 39). 


DISCUSSION 


LAGRANGIAN VARIABLES 


Even a cursory examination of the technique used in this 
investigation for the determination of Lagrangian correla- 
tion shows that, as physical measurements go, this method 
is a “bad” one, largely because of the inherent double 
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(a) \y=0.17. 
(b) £,=0.138. 


FiavurRE 22.—Determination of M and Ly; = 172.3, M. =} inch, and 


U=25.6 feet per second. 


differentiation between measured variable and desired 
information. 

Figures 13 to 22 suggest an uncertainty in values of M 
and L, as large as +20 percent, in spite of moderately good 
precision in the measurement of individual temperature 
distributions such as the lower curve in figure 3. 

As mentioned earlier, the values of »’/U computed from 
initial wake spread are consistently higher than those 
measured with the hot-wire anemometer (reference 18). 


The same relative result was encountered during a brief 


investigation following that reported in reference 18. Up 
to the present time there has been no satisfactory explana- 
tion of the discrepancy. A tentative hypothesis which would 
at least account for its direction may be based upon a human 
weakness in the visual averaging of the reading of a fluctua- 
tion pointer; there seems to be a tendency to choose an 
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9, in. 


(a) Spread of heat from a line source. 
(b) Correlation function R,. 


FiGURE 23.—Spread of heat from a line source and correlation function 


š = 3s 
R, for A{=1 inch, U=8.5 feet per second, and a 38-4- 


“average” more or less halfway between the extremes of the 
needle travel. Thus a pointer motion with very skew prob- 
ability density (greater than 0) would tend to be “averaged” 
at too high a value. The thermocouple voltage in one of 
these thermal] wake traverses has just this character (fig. 11). 
Hence a visual averaging might yield too high a wake width. 
If this effect is nonnegligible, it is advisable to employ some 
electrical means of averaging for skew signals, for example, 
the fluxmeter and bucking circuit described in reference 10. 

In view of the considerable uncertainty in M^ as well as 
that in à and v’ the poor degree of agreement between 
experiment and theory shown in figure 33 is understandable. 


Since the two undetermined constants in the theoretical 


result have been evaluated from sets of experiments com- 
pletely independent of the present ones, this agreement can 
be viewed as an affirmative result. i 

Since some sort of Lagrangian scale should be a significant 
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(a) Spread of heat from a line source. 
(b) Correlation function R,. 
Figure 24.—Spread of heat from a line source and correlation function 


R, for M=1 inch, U=25.6 feet per second, and T5434. 


length in turbulent heat and mass transport, as demonstrated 
in the analytical section of this report, an effort has been 
made to find some systematic variation in the values of L, 
Figure 34 might be construed-to indicate a monotonic 
decrease of L,/L with increasing values of Rz. It is inter- 
esting to note that a decrease was also observed for the 
ratio of mixing length to tube radius by Nikuradse (reference 
30) in fully developed turbulent tube flows. In order to 
determine whether these two rates of decrease with increasing 
Reynolds numbers are of the same order of magnitude, an 
estimate has been made of the magnitudes of Fz correspond- 
ing to Nikuradse’s results given in figures 28 and 29 of 
reference 30. Both scale-to-diameter ratio and average 
turbulent levels for various Reynolds numbers were esti- 
mated with the help of Laufer’s data on turbulent channel 
flow (reference 31) at various Reynolds. numbers. The 
absolute level (i. e., the ordinate scale) of the resulting lmez/ L 
against Rr curve was adjusted to give the most reasonable- 
looking fit with the Z,/Z data. This is the dashed line in 


Ax, it. 


7, in. 


(a) Spread of heat from a line source. 
(b) Correlation function R,. 
Fievurs 25.—Spread of heat from a line source and correlation function 


R, for M=1 inch, Ọ=38 feet per second, and 543.4, 


figure 34 and it shows at least a qualitative resemblance. 
It is likely, however, that £,/L is not a unique function of 
Rr. l 

Of course not all of the scatter in figures 34 and 35 (involv- 
ing L,) can be attributed to simple lack of experimental 
precision. Some is evidence of the fact that the Taylor 
postulate of Lagrangian correlation function being uniquely 
a function of 7 is certainly not very closely true. 
more, table I does show rather systematic variations of L, 
in some of the three-point groups. Most noticeably, there 
is a regular decrease in L, with increasing x,/M (or perhaps 
with decreasing v’/U) for each of the three grids. 


TEMPERATURE-FLUCTUATION FIELD 


Fairly close behind a line heat source in turbulent flow, 
the random pulse nature of the temperature fluctuations at a 
fixed point has been established by the oscillogram in refer- 
ence 10. This is confirmed by the first two oscillograms in 
figure 11, with their highly skew probability densities at 


Further- 
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(a) Spread of heat from a line source. 
(b) Correlation function R,. 


FicuRs 26.—Spread of heat from a line source and correlation function 


R, for M: =7 inch, U=25.6 feet per second, and =9=143.4. 


Azx/\f4=10. One of the objectives of the present investi- 
gation was to find out whether this distinctly pulsed char- 
acter persisted far downstream or whether molecular heat 
conduction becomes increasingly effective in smearing out 
the pulses, until they are no longer distinguishable as such. 
The third oscillogram and probability density in figure 11 
(Ax/4=70) does show a, decided trend away from the pulse- 
type signal. The molecular broadening of the laminar wake 
(corresponding to the pulses) decreases the relative spacing 
of the pulses in g(t) at any point in the “turbulent wake” 
region. ‘This is a reduction in relative length of the flat 
(30) base lines between pulses, giving greater statistical 
symmetry in (t) about its mean, that is, reducing the skow- 
ness of P(®). 

A simple analysis will show the existence of an asymptotic 
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(a) Spread of heat from a line source. 
(b) Correlation function R,. 
FIGURE 27.—S8pread of heat from a line source and correlation function 


Zo —43.4. 


R, for M =- inch, U=25.6 feet per second, and M 


4 


behavior of molecular-conduction effects in a nondecaying 
turbulence. For a nondecaying flow turbulence and very 
large values of ł( =~x/U) the mean-square wake spread due 
to turbulent motion is l 


Y= 2L (79) 


On the other hand, an approach to molecular diffusion 
through Taylor’s concept of “continuous movements” gives, 
for any macroscopic distance downstream, the mean-square 

thermal wake width, 


Ys ASt (80) 


where A is the mean free path and c is the root-mean-square 
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(a) Spread of heat from a line source. 
(b) Correlation function R, 


Fiqaure 28.—Spread of heat from a line source and correlation function 


R, for M=1 inch, V=25.6 feet per second, and =3=86.1. 


molecular velocity. 
From equations (79) and (80) 


i 1/2 
Ta Ac } (81) 


YL v 


For a typical case, take L=1 ‚centimeter, v’=10 centi- 
meters per second, A=6X107" centimeter, and c=5X10* 


7 
m 


centimeters per second. Then Xa =0.17. For people accus- 


tomed to thinking of molecular transport as negligibly small 
in turbulent flow (e. g., in shear flow), this ratio will appear 
quite large. The values of kr/k plotted in figure 38 also 
show that at these low turbulence levels the molecular 
thermal conductivity is not necessarily negligible compared 
with the turbulent transport. 


ʻO yal E re) 4 Re. 
7, in, 


(a) Spread of heat from a line source. 
(b) Correlation function R,. 


FIGURE 29.—Spread of heat from a line source and correlation function 


R, for Af =5 inch, U=25.6 feet per second, and 3g 86.1. 


`~ 


The temperature-fluctuation-level distributions [O across 
the wake (figs. 10 (a) and 10 (b)) show the same character 
as that measured at much higher turbulence level in a jet 
(reference 10), with somewhat lower minimum values, which 
are attributable to the lower turbulence level. A rough 
evaluation of the behavior of the statistical variables in this 
turbulent thermal wake is obtainable by recalling that it 
consists of a randomly ‘waving’ laminar thermal wake. If 
a(t) is crudely represented by a randomly spaced sequence 
of identical rectangular pulses with height #,, width 7, and 
average spacing s, it is easily seen that . 


(82) 


This permits 3/6 to vary between 0 and œ as ¢/j travels the 
permissible range from 1 to œ. Since points nearer the edge 
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(a) Spread of heat from a line source. 
- (b) Correlation function Rẹ 


¥rauRE 30.—Spread of heat from a line source and correlation function 


R, for ar => inch, U=25.6 feet per second, and 57 = 86-1. 


of the turbulent thermal wake have higher values of s/7, the 
behavior of equation (83) is consistent with the experimental 
distribution. Ifthe analysis were repeated with triangular 
pulses, for example, the quantitative estimate would doubt- 
less be more realistic. The higher values of (3//0) min on- 
countered at higher values of v’/U are an indication that 
for a given width of laminar thermal wake, the higher value 
of v’/U leads to a higher minimum value of 8/7. 

The closely Gaussian shape of 6/0, against y has already 
been pointed out. If 


Oe E Er | 
o E aT 
be introduced, there results 
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(a) Spread of heat from a line source. 
(b) Correlation function R, 


Fiaurs 31.—Spread of heat from a line source and correlation function 


R, for M=} inch, U=25.6 feet per second, and 3¢=172.8. 


and 
a e a a Eea yy 68 


Both of these expressions have behavior consistent with the 
experiments. 

The form of dimensionless transverse turbulent-heat- 
transfer rate 30/6,U can be deduced for small values of Az 
(such that R,=1) with this pulse representation of 3). In 
this picture, ga isthe correlation betweeaa continuous Tan: 
dom variable v and a random pulse signal ®© which “fires” 
every time the continuous variable passes through a specific 
value 


Therefore 
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(a) Spread of heat from a line source. 
(b) Correlation function R,. 


FicurB 32.—Spread of heat from a line source and correlation function 


Ra for M= = z inch, U=25.6 feet per second, and a= 172. 3. 


U -= 
oo —— 
ý Az Oy (86) 


where s and © are functions of y. With Gaussian @(y), 


on 
ST ag OP exp | — 2 Trad (87) 


The direct comparison between. this crude picture and the 
experimental results will be confined to the correlation co- 
efficient Rọ = 90/0v'. This is of particular interest in view 
of the surprisingly high experimental values. With equa- 
tion (85) and the fact that »/U=Y’/Az, there results 


Tolman 


This contains the undetermined constant 3/00, which can 
be obtained from any one of several experimental results. 
Figure 36 includes one plot of equation (88) with ¥,/6, de- 


—4 | Ue 
Ras rats exp (8 8) 


--- Equation (35) for A(4)—~O 


“Equation (34) 
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FicurE 33.—Ratio of Lagrangian and Eulerian microscales. 
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FiaurE 34.—Ratio of Lagrangian and Eulerian scales. 


termined from equation (84) and the experimental value of 
(3//O)min and one plot with #,/6 determined by matching 
equation (88) with the experimental result at y/Y’=1.4. 

It is also surprising to find that the experimental (Rss) mer 
at large values of Az is even larger than that at small values 
of Ar. This may be due to a considerable experimental 
error; the resistance-thermometer voltage signal is much 
lower here. Unfortunately no relation corresponding to 
equation (88) has been deduced for large values of Az, 
where R, is essentially zero. 

The criteria. for Taylor’s hypothesis of the interchange- 


ability of instantaneous space U 2 and time 2 derivatives 


(assumed by him to depend only upon turbulence level) 
have been expressed in equations (75) and (78) as functions 
of turbulence level Ry and /A,. If A/d, is replaced by its 
theoretical expression (equation (34)) in terms of R, 
equation (75) becomes 
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Ficure 35.—Ratio of Lagrangian microscale and scale. 
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(b) Re, against y/Y’. Dashed curves are result of analysis based on 
rectangular temperature pulses. 
Ficure 36.—Heat-transfer correlation across thermal wake, computed 


from measured mean temperature distribution. U=25.6 feet per 
second, M=1 inch, z,.=43.4 inches, and Ar=10.5 inches. 
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(a) 6.U against y! 
(b) Re, against y/Y’. 
Fraurs 37.—Heat-transfer correlation across thermal wake, computed 


from measured mean temperature distribution. U=25.6 feet per 
second, M=1 inch, 7.—43.4 inches, and Ar=70 inches. 


v (25 , 29\17 
T VARIA <1 (89) 


For estimates of most flows the first term in the parentheses 
can be neglected; values of 2, less than 5 or 10 are rare. 
Since equation (34) has now been roughly verified by 
experiment, equations (89) and (78) may serve as approxi- 
mate criteria for the validity of Taylor’s hypothesis. 
In the limit of 2,0 when equation (35) replaces equation. 
(34), there follows a simpler criterion to replace equation (75): 


4 [v 

=- | = 1 

3 (F) < (90) 
SUMMARY OF RESULTS 


The following results were obtained from the investiga- 
tion of the diffusion of heat from a line source in isotropic 
turbulence. 
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(a) Turbulent-heat-transfer coefficient at three airspeeds. M=1 inch 
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(b) Turbulent-hent-transfer coefficient for two different grids. 


U=24.6 feet per second and 533 43.4. 


Ficure 38.—Turbulent-heat-transfer coefficient for three airspeeds 
and two different grids, 
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FIGURE 39.—Dimensionless turbulent-heat-transfer coefficient, 
Ze — 43.4 and M=1 inch. 
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1. The thermal wake behind a heated wire set perpendic- 
war to a flowing isotropic turbulence (at sufficiently low 
wire Reynolds number) consists of a randomly “waving” 
thin, laminar, thermal wake whose variations in lateral posi- 
tion give what may be called the turbulent thermal wake. 
At a fixed point not too far behind the wire the instantane- 
ous temperature difference #(é) is a random pulse function, 
and the nature of the turbulent heat transfer can be deduced 
on this basis. Farther downstream the distinct pulse nature 
tends to disappear. _ 


2. The mean transverse temperature distribution = (y) 


appears to be Gaussian within the experimental precision 
for all distances behind the wire. 

3. An Eulerian analysis of this turbulent-heat-transfer 
problem permits computation of the turbulent-heat-transfer 
coefficient Ær which is essentially constant with respect to 
the distance in the direction of the measured diffusion y for 
these boundary conditions. It is found that at low turbu- 
lence levels (approximately equal to 1 to 2 percent) the 
molecular heat transport is not vanishingly small compared 
with the turbulent heat transport. 

4. Although Taylors postulate that Lagrangian correla- 
tions in decaying turbulence can be made similar by intro- 


duction of an independent variable n= Í : v (t)dt (where ¢ is 


time and v’ is the root-mean-square instantaneous veloc- 
ity fluctuation in the y-direction) seems to be an oversim- 
plification, it has been applied here for convenience in the 
reduction of data. A simple comparison of Eulerian and 
Lagrangian analyses for diffusion in nondecaying turbulence 
shows that for large values of the distance from the heat 
source Az the Lagrangian scale L, enters the expression for 
ky, the turbulent-heat-transfer coefficient, much like the 
empirical mixing length in the old turbulent transport the- 
ories. Therefore some properly modified generalization of 
Taylor’s y-postulate should prove useful. 

5. A correction and generalization of Heisenberg’s theo- 
retical expression for the ratio of Eulerian to Lagrangian 
microscale A/A, as a function only of the turbulence Rey- 
nolds number based on microscale Ry has been made and 
seems to agree roughly with experiment. It must be noted 
that since A, depends only upon a transformation dy=v'dt, 
and not upon the integral postulate stated above, its valid- 
ity is not impaired by any failure of the integral postulate 

6. Taylor’s hypothesis for the interchangeability of space 
and time derivatives at low turbulence levels has been ex- 
pressed in terms of criteria which depend upon turbulence 
level, Reynolds number, and A/A,. Applied to the flows 
studied here it shows that in these cases such a transforma- 
tion is permissible. By substitution of the theoretical ex- 


pression for > (Ry), a slightly simpler and rougher criterion 
[r 
is derived, depending only upon turbulence level and Æ. 


TuE Joans HOPKINS UNIVERSITY, 
BALT™ORE, Mb., June 5, 1951. 
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